Eigenmodes of an Asymmetric Cylindrical 

Confocal Laser Resonator with a Single 

Output-Coupling Aperture 

By D. E. McCUMBER 

(Manuscript received February 24, 1969) 

Previous calculations of the low-loss modes of a symmetric cylindrical 
confocal laser resonator have been extended to the asymmetric case. Diffrac- 
tion losses are governed by the geometric-mean Fresnel number N m of the 
two end mirrors and, in the system ice consider, by the Fresnel number N„ 
of an output coupling aperture in one of the mirrors. Loss factors and 
mirror field distributions have been calculated numerically for different N„ 
for N m in the range 0.6 ^ N m ^ 2. 

I. INTRODUCTION 

In a previous paper 1 we described the diffraction losses and the field 
distributions at the reflectors of the low-loss modes of a symmetric 
cylindrical confocal resonator for Fresnel numbers 0.6 ^ N m ^ 2. We 
considered the effect of output-coupling apertures in the reflectors but 
we assumed that both reflectors were identical, each with the same 
output aperture and the same maximum radius. In this paper we again 
consider the cylindrical confocal geometry, but we do not require 
identical reflectors and, in particular, we assume that only one reflector 
is pierced by an output-coupling aperture, as in the coupling scheme 
proposed by Patel and others. 2 

Figure 1 shows an axial section of the confocal resonator in question. 
The cavity is bounded at its two ends by confocal spherical mirrors 
(more exactly, confocal paraboloids 3 ). The first mirror is perfectly 
reflecting over the annular region ^ a l0 ^ p ^ a lm , the second over 
the circular section ^ p ^ a 2m . The maximum radii (a im , a 2m ) are 
both much less than the mirror separation b. 

Expressions for the eigenvalues and eigenfunctions of asymmetric 
rectangular confocal resonators with output coupling slits have been 
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Fig. 1 — Axial section of cylindrical confocal laser cavity. The cavity is 
bounded at its two ends by confocal spherical mirrors with radius of curvature b. 
One mirror is perfectly reflecting over the annular region ^ a i0 ^ p ^ aim, the 
other over the circular section OSp S aa m . Both aim and a 2m are much less 
than b. 

derived by Boyd and Kogelnik. 4 Properties of symmetric resonators 
without coupling apertures are summarized with an extensive list of 
references in the review article by Kogelnik. 5 Equivalence relations 
relating asymmetric and symmetric resonators with circular mirrors 
have been derived by Gordon and Kogelnik. 

Our analysis of the resonator of Fig. 1 closely parallels that of the 
symmetric resonator. 1 Assuming that all dimensions are large com- 
pared with the optical wavelength A, we again use a scalar formula- 
tion of Huygens' principle. 3-5 For the cylindrical confocal geometry 
the field amplitude at reflector ;, j = 1 or 2, for a typical mode can 
be written in the form 

F\i\p, r) = l\i\ P ) exp {-il<p), (1) 

where (p, <p) are radial and angular coordinates in a plane perpendicular 
to the resonator axis and where (I, p) are angular and radial quantum 
integers (transverse quantum numbers). For this asymmetric system 
with nonidentical mirrors, we cannot require for an eigenmode that the 
field amplitude distribution F["(p s <p) at one mirror be a constant 
multiple of that at the other. Rather, we must require for eigenmodes 
that after a round-trip transit of the resonator the field amplitude at one 
mirror be a constant multiple of the initial field amplitude at the same 
mirror. This more elaborate self-reproducing requirement together with 
Huygens' principle gives the following pair of simultaneous integral 
equations which must be satisfied by the radial eigenfunctions j\'J(p) 
and eigenvalues k\^ [compare equation (2) of Ref. 1]. 

"WfiVM = g J dp, p l J l {2Tp 2Pl /b\)j%\p i ), (2a) 
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JfffffW = % [" dp* P2J l (2*p lP2 /b\)f l 2 j :( P2 ). (2b) 

J i (z) is the Bessel function of order \ I \. The loss factor is 

ai v = 1 — | Kijfi" |, (3) 

which is the fractional energy of a mode lost per reflection (or during 
the one-way transit time b/c, where c is the velocity of light in the 
resonator). 5 The phase of the eigenvalue product k$k$ determines the 
resonant wavelength: 

resonant X - 4tt6/[(Z + 1)tt - Arg icgitf? - 2tto], (4) 

where n is an arbitrary integer (longitudinal quantum number) . 
It is useful to introduce the mirror Fresnel numbers 

tf™ = a\J\b, Ni 2) = a\J\b, (5a) 

and their geometric mean 

N m e- [N™N™]* = a lm a 2m /\b. (5b) 

In place of the variables p, and the functions f£(j> { ) in equations (2), 
we introduce new variables 

fi = r m pj/a im (6a) 

and functions 

g\i\r 4 ) = tii\r iai Jr m ), (6b) 

where r 2 = N m . We characterize the size of the radius-a,,, hole in 
the first mirror by the Fresnel number 

N m r» = {aJa lm ) 2 N m = a 2 a 2m /\ba lm . (7) 

With no significant loss of generality we can define the functions /{J } (p,) 
such that 

a 2m K\V/a lm = a lm K\l ) /a am = k 1v (8) 



and 



5 pa = 2tt J " dr x r^ir^giHh), (9a) 

«„ = 2r f dr>r 2 gZ\r 2 )g? a \r 2 ). (9b) 

» a 
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With this notation the eigenvalue equations (2) become 

«cippi?(r a ) = 2tt f" dr, r 1 S l Chf t r l )£'fa)i (10a) 

Ki^S'fri) = 2t f dr 2 r 2 J l (2irr 1 r 2 )g\ 2 p ) (r a ) . (10b) 

Equation (3) simplifies to 

Clip = 1 — | *lv \ 2 > (11) 

which depends upon the parameters (a l0 , a lm , a 2m , b, X) only through 
the Fresnel numbers N m = r* and N = r\ . In what follows we de- 
scribe how the loss factor a lp and the amplitudes g\ p {r) change with 
N for N m in the interval 0.6 ^ N„ ^ 2. Solutions for N = are 
described elsewhere. 1 ,3,7 ~ 9 

Our numerical method is similar to that previously described for 
the symmetric resonator. 1 We expand the Bessel-function kernels 
in equation (10) as power series, truncate the series after a finite 
number M = max (10 N m + 1, 10) of terms, and reduce the integral 
equations (10) to M-dimensional matrix equations which are solved 
numerically with standard matrix routines. [The reduction of equa- 
tion (10) to matrix form is described in the Appendix.] The merits 
of this technique remain as described before. 1 

II. ANALYTIC METHODS FOR SMALL JV„ 

The eigenvalues k, p and field amplitudes g tp (r) for N e = are de- 
scribed elsewhere. 1 ' 3 7_0 The loss factor a lv = 1 — | k 1p | 2 for the four 
lowest-loss modes (compare with Fig. 2 of Ref. 1) are tabulated for 
0.6 ^ N m ^ 2 in Table I. For I = 0, the field amplitude g lv (r) at r = 
is finite; for I 7*- 0, it vanishes as r ' " . As before, we expect the modes with 
angular quantum number I = to be more sensitive to the coupling 
aperture (N„ > 0) than the I ?* modes. 1 

If we use a superscript "0" to lable the eigenvalues and field ampli- 
tudes for N„ = 0, then for small r, to within corrections of relative 
order r 2 , 

9°M = g° 0p (0) for 1 = (12a) 

= c? p r'" for I 9*0, (12b) 

where coefficients gl P (0) are listed in Table II (an expanded version of 
Table III, Ref. 1) and coefficients c% in Table III. The forms (12) are 
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Table I — Loss Factor a, p for N — 



Nm 


aoo° 


aoi° 


aio° 


OI0° 


0.6 


3.614 


X 


10" 1 


0.7931 




0.2724 




0.6679 


0.7 


1.301 


X 


10-* 


0.6131 




0.1371 




0.4712 


0.8 


4.448 


X 


io-» 


0.4131 




6.103 X 


io-» 


0.2900 


0.9 


1.471 


X 


io-» 


0.2411 




2.477 X 


io-« 


0.1563 


1.0 


4.759 


X 


10-* 


0.1233 




9.417 X 


IO" 3 


7.505 X IO" 2 


1.1 


1.515 


X 


io-« 


5.651 X 


io-» 


3.424 X 


io-» 


3.285 X IO" 1 


1.2 


4.767 


X 


10" 6 


2.382 X 


io-» 


1.206 X 


IO" 3 


1.343 X IO" 2 


1.3 


1.485 


X 


io-« 


9.462 X 


10~ J 


4.151 X 


10"* 


5.225 X IO" 3 


1.4 


4.59 


X 


io-« 


3.601 X 


10"» 


1.403 X 


io-* 


1.961 X IO" 3 


1.5 


1.41 


X 


io-« 


1.328 X 


lO- 3 


4.674 X 


10-* 


7.164 X 10-* 


1.6 


4.3 


X 


10~ 7 


4.78 X 


io-* 


1.539 X 


10~ B 


2.561 X 10"* 


1.7 


1.3 


X 


10" 7 


1.69 X 


io-* 


5.01 X 


IO" 6 


9.000 X IO" 5 


1.8 


4 


X 


io-» 


5.89 X 


10" 6 


1.62 X 


IO" 6 


3.116 X IO" 5 


1.9 


1 


X 


io-» 


2.02 X 


io- B 


5.2 X 


10~ 7 


1.066 X IO" 6 


2.0 


4 


X 


io-» 


6.86 X 


io-« 


1.7 X 


IO"' 


3.60 X IO* 8 



useful for estimating the perturbations induced by a small finite N„ . 
To first order, the perturbed field amplitudes are 



~ E' g°M /j. J K|,) /ov» 2tt f" dr x wlMglM, (13a) 

Q*V \ K lv) \ K lQl J 

Table II — Field Amplitude at r = for I = 
Modes with N B = 



Nm 


&00°(0) 


eo.°(0) 


(70>°(O) 


0.6 
0.7 
0.8 
0.9 


1.2770 
1.3021 
1.3213 
1.3354 


1.2251 
1 . 1541 
1 . 1254 
1 . 1286 


1.6159 
1.4851 
1.3735 
1.2760 


1.0 
1.1 
1.2 
1.3 


1.3457 
1.3536 
1.3597 
1.3647 


1.1511 
1 . 1807 
1.2093 
1.2336 


1 . 1930 
1.1291 
1.0890 
1.0742 


1.4 

1.5 
1.6 

1.7 


1.3688 
1.3723 
1.3752 
1.3778 


1.2532 
1.2688 
1.2814 
1.2918 


1.0812 
1 . 1025 
1 . 1302 
1 . 1580 


1.8 
1.9 
2.0 

CD 


1.3800 

1.3820 

1.3838 

■v/2 = 1.4142 


1.3006 
1.3081 
1.3146 
1.4142 


1.1826 
1.2033 
1.2203 
1.4142 
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Table III — Field-Amplitude Coefficient for 
I -^ Modes at r = for N = 



N m 


C10° 


cm" 


C20° 


0.6 
0.7 
0.8 
0.9 


2.6428 
2.7222 
2.8269 
2.9266 


4.4649 
4.0086 
3.7024 
3.5247 


3.9948 
3.9131 
3.9979 
4.1770 


1.0 
1.1 
1.2 
1.3 


3.0094 
3.0746 
3.1255 
3 . 1659 


3.4618 
3.4932 
3.5889 
3.7148 


4.3917 
4.5999 
4.7809 
4.9305 


1.4 
1.5 
1.6 

1.7 


3 . 1988 
3.2260 
3.2491 
3.2690 


3.8431 
3.9580 
4.0549 
4.1351 


5.0526 
5.1531 
5.2371 
5.3087 


1.8 
1.9 
2.0 

CO 


3.2862 

3.3014 

3.3148 

2w in = 3.5449 


4.2017 

4.2578 

4.3059 

23/vi/a = 5.0133 


5.3704 

5.4244 

5.4721 

2tt = 6.2832 







•2ir f" drxTxfudrdfM- ( 13b ) 

To second order, the loss factor is 



«i» 



- a° lp + (1 - a° lp )2T f° dr x fJMrdf 

•'o 

_ E , (1 ~ a?.)(l - eel) L f' ^ irig o Mg o M li m 
fS. a lq - a lv L Jo J 



(14) 



In deriving (14), we used the fact that for the cylindrical confocal 
geometry the eigenvalues k [p are real and k\ v = \ k 1v \ . 

The last terms in (13) and (14) describe mode mixing by the aper- 
ture. The amount of mixing depends upon the separation of the 
eigenvalues as well as upon the strength of the perturbation. Degener- 
ate or nearly degenerate modes are much more sensitively coupled than 
are modes with greatly different losses. In the symmetric resonator the 
even-p and odd-p modes of a particular angular quantum number I 
do not mix; such modes do mix in the asymmetric geometry. 1 

For N sufficiently small, we can neglect the second-order or mode- 
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mixing term in (14). If we use the small-r approximations (12) in the 
remaining integral, we find for I = that 

* 0P = «o° P + (1 - *o°>tfo[<&(0)] 2 (15a) 

and for 1^0 that 

«, p = al + (1 - a ?,M4W" M /(| « I + I)- (15b) 

These expressions confirm our previous conjecture that modes with 
I = are more sensitive to aperture loss than are modes with I ^ 0. 
A quantity of interest in the design of lasers with aperture output 
coupling is that value N oc of N for which the losses of the longitudinal 
(00) mode equal the losses of the (least lossy) transverse (10) mode. 2 
All other things being equal, the laser will operate in the (00) mode 
for iV < N oc , in the (10) mode for N > N oc , and in still another 
mode for larger values of N . Using Eqs. (15), we estimate that 

n bc = («;„ - ao° )A(i - <4«o(o)] 2 . (i6) 

III. NUMERICAL RESULTS FOR N FINITE 

Let us compare estimates based upon the approximate expressions 
(15) and (16) with accurate numerical results. Using the numerical 
technique outlined at the end of the introduction and in greater detail 
in the Appendix, we have computed the loss factor a /p for N finite and 
N m in the range 0.6 ^ N m ^ 2. Results for N m = 0.8 and N variable 
are shown in Fig. 2 and similar results for iV m = 1.6 in Fig. 3. Results 
for N„ = 0.001 and N m variable are shown in Fig. 4. These examples 
were chosen to facilitate comparison with the two-aperture symmetric 
geometry of Ref. 1. The results are qualitatively similar to those ob- 
tained before, except that here odd-p and even-p modes interact whereas 
before they did not. 

Predictions based upon the first-order expressions (15) and Tables 
I to III are shown as dashed lines in Figs. 2 and 3. The fit to the exact 
results is good for sufficiently small N , but deviations are large for 
those A r 's for which the interaction between the {I, p) and (I, p + 1) 
modes is evident as a repulsion in the calculated loss curves. 

The critical single-aperture Fresnel number N oc for which the loss 
factor ai p of the longitudinal (00) mode equals that of the lowest 
transverse (10) mode is shown as a function of mirror Fresnel number 
N m in Fig. 5. This is an important parameter in the design of aperture- 
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Fig. 2 — Loss factor ai P versus single-aperture Fresnel number N for low-loss 
modes when N m = 0.8. (Compare with Fig. 9 of Ref. 1 for the symmetric two- 
aperture geometry.) The dashed lines are estimates based on equation (15) and 
the data from Tables I through III. 

out-put-coupled cavities having good mode selection. 2 Also shown in 
Fig. 5 is the estimate of N oc derived from Eq. (16) and the data from 
Tables I through III. The agreement is reasonably good, and consistent 
with what one would expect from the accuracy of the estimates derived 
from Eqs. (15) in Figs. 2 and 3. 

The effect of mode coupling is apparent in the amplitude g\p(r) of 
the field at the two mirrors. Consider the 1 = modes, which are those 
most sensitive to finite N„ . The field amplitudes and intensities at the 
mirrors of the two lowest-loss modes are shown for N m = 0.8 in Figs. 6 
and 7 and for N m = 1.6 in Figs. 8 through 11. Figures 6 and 8 show 
the distributions for N = 0; they are the same on both mirrors. [Distri- 
butions for other (Ip) modes are shown for N = in Ref. 1.] The other 
figures show how the distributions change for N > 0. In each figure 
the dashed lines indicate the field distributions on the mirror pierced 
by the aperture, the solid lines those on the intact mirror. The radius of 
the mirrors and the radius of the aperature are indicated on the plots. 

One should distinguish two effects apparent in the field plots as N B 
increases. First, there is a change in the magnitude of the field ampli- 
tude g™(r) on the pierced mirror. This is a simple renormalization 
correction implicit in the requirement (9) that the fields be normalized 
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over the reflecting areas of the mirrors. Second, there are changes in 
the shape of the field distributions on both mirrors. This is a consequence 
of mode mixing, which becomes appreciable for those N 's for which 
significant mode repulsion is apparent in the curves of Figs. 2 and 3. 
In each case the effect is to reduce the intensity of the less-lossy mode 
in regions where the mirrors are not reflecting, at the expense of the 
more lossy of the two interacting modes. This is apparent, for example, 
in Fig. 9 with N m = 1.6 and N = 0.0001, for which there is appreciable 
interaction between the (00) and (01) modes (compare with Fig. 3). 
The amplitude of the (00) mode at the aperture is decreased below that 
in Fig. 9; that of the (01) mode is increased. In Fig. 11 with N m = 1.6 




10~ J 10" 

APERTURE FRESNEL NUMBER Nj' 2 ' 



Fig. 3 — Loss factor ai p versus single-aperture Frcsnel number A r „ for low-loss 
modes when JV m := 1.6. (Compare with Fig. 11 of Ref. 1 for the symmetric 
two-aperture geometry.) The dashed lines are estimates based on equation (15) 
and the data from Tables I through III. 
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0.6 0.8 1.0 1.2 1.4 1.6 I. 

APERTURE FRESNEL NUMBER N m 

Fig. 4 — Loss factor a« P versus Fresnel number iV m for low-loss modes with No 
fixed at 0.001. (Compare with Fig. 16 of Ref. 1 for the symmetric two-aperture 
geometry.) 

and N = 0.01, the amplitude of the (01) mode at the aperture is re- 
duced as a consequence of interaction with the (02) mode (Fig. 3). 



IV. DISCUSSION 

Our previous treatment 1 of a symmetric cylindrical confocal laser 
cavity is extended here to the asymmetric cylindrical confocal geometry 
of Fig. 1 and specific numerical values for the loss factor and for the 
field distributions at the mirrors have been calculated for Fresnel 
numbers 0.6 ^ N m ^ 2. The transformations outlined in Section I 
show that the loss factor of a cavity with different-sized mirrors 
(N^ 9^ N™) equals that of a cavity having two mirrors with the 
same outer dimensions [N m = (N^N™)*]. The field distributions 
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scale accordingly [see equations (6)]. Similar results also obtain in the 
rectangular confocal geometry. 4 

For sufficiently small aperture Fresnel numbers iV the cavity losses 
associated with a single output coupling aperture can be estimated as 
in equations (15) from first-order perturbation theory. Just as in the 
symmetric two-aperture case considered previously, the value of N 
for which such first-order calculations fail decreases rapidly as the 
Fresnel number N m increases, because the field distributions distort 
through mode mixing to minimize the aperture losses in the lowest- 
loss modes. 1 As before, this distortion occurs at approximately those 
values of N and iV m for which an observer at one reflector, using light 
of the relevant wavelength and optics limited by the radius r m = N£ , 
can resolve the aperture of radius r = N\ from the other reflector. 1,10 

APPENDIX 

Reduction of Integral Equations to Matrix Equations 

We express the Bessel-f unction kernels in power-series form. (I = 
I 1 1 throughout this appendix.) 



\\ 


\\ 
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\w EXACT 
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/ \ 


y \ 


ESTIMATED \ 
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...\. 



0.8 1.0 1.2 

FRESNEL NUMBER Nr 



Fig. 5 — Critical single-aperture Fresnel number A^oo for which diffraction 
losses of longitudinal (00) mode equal those of the lowest transverse (10) mode 
versus Fresnel number N m . (Compare with Fig. 18 of Ref. 1 for the symmetric 
two-aperture geometry.) The dashed line is an estimate based on equation (16) 
and the data from Tables I through III. 
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Fig. 6— (a) Field amplitude g tp (r) and (b) field intensity | g, P (r) | 2 for modes 
dp) = (00) and (01) with N m = 0.8 and N = 0. The field distributions are 
identical on both mirrors. 



J z (27rr 1 r 2 ) = (7rr,r 2 )' X) 



-1/ „ - \2(m-l) 



(-l)"- 1 ^,^) 



%(to + I- 1)!(to- 1)! 



(17) 



Truncating this series after M terms and substituting the result into 
(10a), we obtain 



dr l r l X) 



c-ir-w a ) 



l+2(m-l) 



(to + Z - 1)! (to - 1)! 



fiffto (18a) 



[(to- 1)!(to+ I - 1) !//!]* m KhP) ' 



(18b) 



where 
G2\l, v) = 



2tt 



[(to + Z - 1)! (to - 1) 



^ [ n dr 1 r l {rrT-^ ln g\l\r,). 



(19) 
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Likewise we obtain from (10b) 

«^?w = pfO' t [(m J;Z m Yf-wiy G '"«' *■ (20) 

where 

Substituting the expression (20) for g\" (r t ) into the right hand side of 
(19) and the expression (18b) for gfffa) into (21), we obtain after 
simple manipulations 

k 1p G1 1 \1, p) 

iT^w'T,: wrn , , n ^(i, P)l 



., [(m-1)! (m+Z-1)! (fc-1)! (fc+Z-l)!]*(H-ro+A:-l) 



(22a) 




Fig. 7— (a) Field amplitude g tp <n (r) and (b) field intensity | pi P (i> (r) | 2 for 
modes (Zp) = (00) and (01) with N m = 0.8 and N. = 0.03. The dashed fines 
refer to mirror 1 (Fig. 1) and the solid lines to mirror 2. The radius of the 
aperture in mirror 1 is r = 0.173. 
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RADIUS r 



Fig. 8— (a) Field amplitude g, P (r) and (b) field intensity | g lp (r) | 2 for modes 
dp) = (00) and (01) with N m = 1.6 and N* = 0. The field distributions are 
identical on both mirrors. 



K lp Gl 2) (l,p) 



= z 



(-ir\TN m y +m 



l+m+k-l 



\ [(m-1)! (m+7-1)! (fc-1)! (fc+/-l)!]*(H-m+A;-l) 



Gi l) (l, P ). 

(22b) 



We have used the definitions (5b) and (7) to replace (r* , r 2 ) by the 
Fresnel numbers (N m , N ). 

It is convenient to view G" y (l, p) as the mth component of an M- 
dimensional vector G U) (l, p). We define B to be the M X M diagonal 
matrix with elements B mm = (— l)"" 1 . We also define real symmetric 
matrices S ( "(Z) and S (2) (Z) with elements 



OQ = 



(wN m ) 



l+m+k-l 



[(m-1)! (m+Z-1)! (fc-1)! (fc+Z-l)!]*(H-m+ft-l) 

[(*N m y +m+h - 1 - fay.)'*"**" 1 ] 

[(m-1)! (m+Z-1)! (fc-1)! (fc+«-l)!]*(/+m+fc-l) 



, (23a) 
■ (23b) 
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With these definitions equations (22) can be written more compactly 



k„G (1) (Z,p) = S (a) (Z).B.G (2 Up), 
K lp G l2, (l,p) = S ll \Q>B'G ll) (l,p). 
Eliminating G (2> , we obtain 

KlG ll \l,p) = S (2, (/)-B-S u) (0-B.G (1 Up), 



(24a) 
(24b) 

(25) 



which is a single il/-dimensional matrix eigenvalue equation. 

It is generally useful to transform equation (25) such that the matrix 
to be diagonalized is real symmetric (or Hermitian). Because S ( ~'(/) 
is real symmetric with nonnegative eigenvalues, we can find a real 
lower-triangular matrix P(/) such that 



S (2) (/) = P(/)-P(/) r . 



(26) 




Fig. 9 — (a) Field amplitude ffi P 0) (r) and (b) field intensity \gi P u> (r) \- for 
modes (Ip) = (00) and (01) with A r m = 1.6 and N = 0.0001. The dashed lines 
refer to mirror 1 (Fig. 1) and the solid lines to mirror 2. The radius of the 
aperture in mirror 1 is r„ = 0.01. 
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If we define a new vector 

F(l,p) = P(0" 1 -G U) a P ), 



(27) 



then (25) can be written 

K 2 IP V(1,P) = P(Z) r -B-S (1) (0-B-P(0-F, (28) 

for which the matrix on the right hand side is obviously real symmetric. 
If U(0 is the real orthogonal matrix which diagonalizes this matrix, then 

U(Q-K(Z) = P(Z) r -B-S (1) (Z)-B-P(Z)-U(Q (29) 

where K(l) is diagonal with elements K„(J) = k] p , V — 1 to M. The 
eigenvector F(Z, p) of (28) corresponds to the pth column of V(l) and, 
from (27), 

O m (l,p) = P(Z)-F(Z,p). (30) 




Fig. 10— (a) Field amplitude g tp w (r) and (b) field intensity \g, r l,i (r) | 2 for 
modes Up) = (00) and (01) with N m = 1.6 and N, = 0.001. The dashed lines 
refer to mirror 1 (Fig. 1) and the solid lines to mirror 2. The radius of the 
aperture in mirror 1 is r. = 0.0316. 
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Fig. 11 — (a) Field amplitude g lp w (r) and (b) field intensity \gi, w (r) | 2 for 
modes (Ip) = (00) and (01) with N m = 1.6 and No = 0.01. The dashed lines 
refer to mirror 1 (Fig. 1) and the solid lines to mirror 2. The radius of the 
aperture in mirror 1 is r = 0.1. 

The elements of U(Z) and K(l) are easily computed numerically; the 
vectors G (1) (Z, p) follow from (30); the vectors G (2) (Z, p) follow from 
(24b); and the amplitudes g\l\r) follow from (18b) and (20). 

The program used to compute the results reported in this paper 
required a nominal 0.0003 hr. of GE 635 processor time to compute the 
M different eigenvalues | k 1v | and eigenvectors [G (1) (Z, p), G (2) (l, p)] 
for M = 20. Timing for other values of M varies roughly as M 3 . 
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